arXiv:1509.06206vl [cond-mat.supr-con] 21 Sep 2015 


Magnetization of underdoped YBa 2 Cu 30 y above the 

irreversibility field 

Jing Fei Yu/’Q B. J. Ramshaw,^ 1. Kokanovic,^’^ K. A. Modic,^ 

N. Harrison,^ James Day,® Ruixing Liang,®’® W. N. Hardy,®’® D. 

A. Bonn,®’® A. McCollam,^ S. R. Julian,^’® and J. R. Cooper® 

^Department of Physies, University of Toronto, 

Toronto, Ontario, MBS 1A7, Canada 
‘^Los Alamos National Laboratory, Mail Stop E536, 

Los Alamos, New Mexico, 815f5, USA 
^Cavendish Laboratory, University of Cambridge, 

Cambridge, CBS OHE, United Kingdom 
‘^Department of Physics, Paculty of Science, 

University of Zagreb, P. O. Box 331, Zagreb, Croatia 
^Department of Physics and Astronomy, 

University of British Columbia, Vaneouver, 

British Columbia, V6T IZl, Canada 
^Canadian Institute for Advanced Research, 

180 Dundas St.W, Toronto, Cntario, MBS 1Z8, Canada 
High Pield Magnet Laboratory, Radboud University, 
Toernooiveld 1, 6B2B ED Nijmegen, the Netherlands 
(Dated: August 18, 2016) 


1 


Abstract 

Torque magnetization measurements on YBa 2 Cu 30 y (YBCO) at doping y = 6.67(p = 0.12), in 
DC fields (B) up to 33 T and temperatures down to 4.5 K, show that weak diamagnetism persists 
above the extrapolated irreversibility field = 0) ~ 24 T. The differential susceptibility 

dM/dB, however, is more rapidly suppressed for B > 16 T than expected from the properties of 
the low field superconducting state, and saturates at a low value for fields B > 24 T. In addition, 
torque measurements on a p = 0.11 YBCO crystal in pulsed field up to 65 T and temperatures 
down to 8 K show similar behaviour, with no additional features at higher fields. We discuss 
several candidate scenarios to explain these observations; (a) superconductivity survives but is 
heavily suppressed at high field by competition with CDW order; (b) static superconductivity 
disappears near 24 T and is followed by a region of fluctuating superconductivity, which causes 
dM/dB to saturate at high field; (c) the stronger 3D ordered CDW that sets in above 15 T may 
suppress the normal state spin susceptibility sufficiently to give an apparent diamagnetism of the 
magnitude observed. 

PACS numbers: 74.72.Gh, 74.25.Ha, 74.25.Op, 74.25.Bt 
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The possible existence of bulk superconductivity as T —)• 0 K above the irreversibility 


held (-ffirr)lll in the cuprates has been a long standing question. Not only is this problem 
important for our understanding of the cuprates, but also because there is still debate}^, 
Q] about whether Cooper pairs persist in the region of the held-temperature plane where 
quantum oscillations are seenjd]. 

Many experimental efforts have been made to address this issue [ 
has consistently been reported using torque magnetometry at high helds in many families 
of cuprates and it is argued that this observation shows the persistence of Cooper pairs 


Diamagnetism 


above 


irrfl. 


For YBa 2 Cu 30 j^, resistivity measurements have established = 0) to 


be below 30 T for helds along the c-axis for c 


(OVIII)Q. 


Moreover, X-ray 


10 


opings between p = 0.11 (Oil) and p = 0.12 
, lll|, NMR 12|, and sound velocity measurements [l3| have 


demonstrated the existence of static charge density wave (CDW) order that competes with 
superconductivity: Ref. [m| shows a distinct long range 3D order that continues to grow at 
28 T for an OVIII crystal. The CDW is strongest and the suppression of Hc 2 is largest at 
p = 0.125 for YBCO[ld, [ij. 

Recent thermal conductivity measurements by Grissonnanche et al. show a sharp 
transition precisely at the extrapolated Trirr(T = 0) ~ 22 T for Oil YBCO. They have 
interpreted this feature (henceforth referred to as Hk) as a signature of Hc 2 , arguing that 
the end of the rapid rise in thermal conductivity at 22 T rehects a corresponding increase 
in the mean-free-path as a result of the sudden disappearance of vortices at Hc 2 - On a 


crystal with same doping, Marcenat et al. [15[ show that the specific heat saturates at a 
field Hep. Hcp(T) lies above but they extrapolate to the same value at T = 0 Kjl5|. 

In contrast, torque measurements by F. Yu et al. jo] on a crystal with the same doping 
suggested diamagnetism persisting to helds much higher than Hk- The debate is thus still 
open. 

To resolve this problem, we conducted torque magnetometry measurements of magneti¬ 
zation (M) on two p = 0.12 (OVIII, = 65 K) crystals in DC helds and one p = 0.11 (Oil, 
Tc = 60 K) crystal in pulsed helds. The crystals were mounted on piezoresistive cantilevers 
and placed on a rotating platform, with the Cu 02 planes parallel to the surface of the lever. 
DC held sweeps, hrst from 0 to 10 T and later from 0 to 33 T, were performed with the c-axis 
of the OVIII crystal at a small angle 6 from the held. The magnetoresistance of the levers 
was eliminated by subtracting data from the complementary angle {—0). Similar procedures 
were used for the Oil crystal in pulsed magnetic helds up to 65 T. For strongly anisotropic 
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FIG. 1. (Color online) Black dots: high temperature anisotropic susceptibility xd{T) of the OVIII 
crystal at 10 T. Blue solid line: fit to this data above 120 K using Eq. [TJ The parameters A = 11.09 

16], while the fit gives = 5.84 A/m/T , T* = 330 


A/m/T, Tp = 680 K are taken from Ref. 


K and X'^(O) = 1-26 A/m/T; Red dashed line: Linear fit with xi^) = 1-22 x 10 ^ x (T + 948), 
following Ref. ^ but with different parameters. Note that 1 x 10“^ emu/mol= 9.73 A/m/T. 

superconductors, where out-of-plane screening currents can be neglected, the torque r per 
unit volume V at an angle 6 from held B is given by r jV = ^xd{T)B‘^ sin 29 +McB sin^jl?]. 
Here Xd(T) = Xc(T) — Xab(T) is the anisotropic susceptibility in the normal state and Me 
is the magnetization from in-plane screening currents for a held of B cos 6 along the c-axis. 
This is a good approximation when Me 3> XdB or when the superconducting gap and Me 
are both small. 

A key challenge with magnetization measurements in the cuprates is the separation of the 
normal state from the superconducting contributions, because superconducting huctuations 
are thought to contribute to x(T) even at temperatures far above Tc 17|, while jg 

temperature dependent to well below Te- We follow the procedure outlined in Refs.[^, 17] 
and interpret Xd{T) in the normal state of underdoped YBCO as arising from the pseudogap 
and ^-factor anisotropy, plus a superconducting huctuation term that sets in below 120 K. 


Neglecting isotropic Curie and core suscepti 


total normal state contribution to Xd(T) is m-- 


Dility terms, which do not contribute to r, the 


x”'(T) = xg°(T) + x'd' + xg(r) 


( 1 ) 


where xY)^ is the T-independent Van Vleck susceptibility Xn^iT) is the pseudogap con¬ 
tribution assuming a V-shaped density of states (DOS)[l8|, and Xoi'^) is thought to arise 
from an electron pocket or Fermi arcs in the region 0.0184 < p < 0.135. Specifically, 
Xd^ = a (1 — y~^ In [cosh( 2 /)]), where A = Vo/r|, y = Eg/2kBT, Eg = ksT* and T* is the 
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pseudogap temperature, and Xd(^) = X'^(O) [1 “ exp(—Tp/T)] where Tp is the Fermi tem¬ 
perature. The £t is shown in Fig. [H along with a linear model for the normal state x used in 
Ref. jo]. Both fits agree well with the data for T > 120 K. Our background is almost twice 
as small as that of the linear fit at T = 0 K. Subtraction of the background magnetization 
using this non-linear model should thus give a signihcantly weaker diamagnetic signal at 
T —)■ 0 K than the linear model would (about 160 A/m at 33 T). 


In Figs. 2(a) and 3(a), we show vs curves at selected temperatures for the OVIII 
and on crystals, obtained by subtracting Mbg = XbgB, where xbg is fhe blue line in Fig. 
m and Bz is the held projected along the crystalline c-axis. For the OVIII crystal, at T = 103 
K, we see that Me is almost zero. At 58 K, just below Tc, we see significant diamagnetism 


that gradually tends to about —130 A/m at high field. Fig. 2(a) shows that the crystal 
remains weakly diamagnetic down to 4.5 K in fields up to 33 T. Similar behaviour was found 


for the on crystal in pulsed helds. As shown in Fig. 3(a), is still diamagnetic at the 
highest held B^ = 63 T, but has a small value - about —90 A/m at 8 K. Our results differ 
from those of F. Yu et al. jo]: for example, at 10 K and 20 T we measure Me to be four 


times larger [19|]; furthermore, at 30 T we hnd about -200 A/m for Oil and OVIII rather 
than their value of -75 A/m. Our estimated uncertainty in Xp)(0) corresponds to ±32 A/m 
in Me at 33 T and ±61 A/m at 63 T. 

Although the weak diamagnetic signal persists to higher helds, we are able to see a 
signature in our diherential susceptibility dM/dB at helds comparable to Hk (22 T) found 
by thermal conductivity j^. In each curve of Fig. 2(b) and 3(b), dM/dB decreases linearly, 
up to a held we call Hd{T), before saturating to a small but non-zero value. At the lowest 
temperatures for both OVIII and Oil crystals, we hnd Ha ~ 24 T, which is close to the 
extrapolated Hi„{T = 0). This is consistent with the feature at found by thermal 
conductivity!?], though unlike Mp, Hd does not correspond to a sharp transition. Hd varies 
very little with temperature for T < 10 K, a result that is consistent with the hndings of 


Ref. 


0 . 


though the T-dependence at high temperatures is not consistent with that found 
by Refs. |^, 15|. Surprisingly, we do not observe in any of our crystals the broad peak in 
dM/dR reported by Ref. [6]. 

In highly anisotropic type-II superconductors, the magnetization calculated using mean 
held (MF) Ginzburg-Landau (GL) theory for an s-wave superconductor, which we use in 
the absence of a d-wave theory, yields logarithmic behaviour at low held (in cgs units), 
—AttM = a(/)o/(87rA^) ln(/9idc2/M) for 0.02 < H/He 2 < 0.3, where a and [3 are numbers of 
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FIG. 2. (Color online) (a) Magnetization (Me) of the OVIII crystal vs B^, the held parallel to the 
c-axis. Here Mc{T, H) = M^hsiT^ H) — Mbg{T, H), where Mbg = XdB and xd is the blue line in 
Fig. [TJ Dashed line: Mbg at 4.5 K. Diamagnetism is present even at our highest held of 33 T. |(b)| 
Differential susceptibility dM/dB of the OVIII crystal vs B^ at selected temperatures. The lines 
are guides to the eye. We call the characteristic held at which dM/dB departs from linearity Hd- 
Red: calculated mean held dM/dB near Hc 2 with k = 50, with k = 41 (purple) and with k = 150 
(blue). 


order 1, 0o is the flux quantum for Cooper pairs and A is the London penetration depth 201. 
/iSR at low fields has shown a ^/H field dependence for A(T = 0) 2l|, but results of tunnelling 


experiments on Bi-2212 imply thermally induced pair breaking near the nodes 


22j, indicating 


a weaker field dependence at higher T. Thus, for simplicity, we assume a negligible held 
dependence of A. We also assume a = 0.77 and /d = 1.44 for 0.02 < H/Hc 2 < 0.3 in 
reasonable agreement with later works 23|, [2^, and we ht the low held magnetization and 
obtain an estimate of Hc 2 (T), shown in Fig. [5l Since our GL values of Hc 2 join smoothly to 
Hd, it is possible to interpret Hd as the low temperature GL type Hc 2 . 
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(a) 



(b) 


FIG. 3. (Color online) (a) Magnetization (Me) of the Oil crystal measured in pulsed magnetic 
field up to Bz = 63 T, where Me = Mobs — Mbg-, Mbg = XdB and xd is the blue line in Fig. 
[TJ For clarity only the falling-field sweeps are shown. Diamagnetism is present though extremely 
weak at high field (inset). The small offset in Me between the T < 40 K and T > 50 K curves may 
be due to the transition to l ong -range CDW order near 40 K in high fields as observed in both 
sound velocity [l^ and NMR 12|. |(b)| Differential susceptibility for the Oil crystal in pulsed field. 
dM/dB is seen to be small and constant up to the highest field of 63 T. Blue: calculated mean 
field dM/dB near He 2 with k = 50. 


When H/He 2 > 0.3, and again using cgs units for an s-wave superconductor, the magne¬ 
tization is expected to obey dvrM = {H — / (2k^ — 1)/?^], where k is the GL parameter 


and Pa = 1-16 is the Abrikosov parameter 2^, l26j. Figs. 2(b) and 3(b) show that for B > 28 
T, dM/dB has the mean field property of saturating toward a constant value, but this is 


very small and requires k ~ 150, a value far greater than k = 50 given Ref. 


3 - 


This means 


that our high held dM/dB is nearly ten times smaller than would be expected. This may 
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FIG. 4. Magnetization data of the OVIII crystal at 16 K. The blue dashed line shows the MF 
behaviour near Hc 2 for an s-wave superconductor with k = 41. The stronger (3D) CDW sets in 
above 15 T for OVIII YBCO. At higher fields the data are consistent with n = 145(soZi(i line). 


be due to the field dependent charge density wave (CDW) order within the vortex liquid 
region^, 11]. The CDW competes with superconductivity and is partially suppressed at 
low field. As increasing field suppresses superconductivity, the CDW order is gradually 
restored ^ . The presence of a relatively strong CDW would increase A and thus increase k, 


as illustrated in Fig. 01 A linear region in Mc{B) can also be seen in Fig. 2(a), for T = 20 


K and T = 16 K and i? < 17 T, with k, = 41, and in Fig. 3(a), for T = 20 K and B < 17 
T, with K = 50. These linear regions are not present above 20 K, where Mc{B) is likely to 
be smeared out by thermal fluctuations. As shown in Fig. 01 for the OVIII crystal, low field 
Me extrapolates to zero around 24 T, consistent with our GL type Hc 2 . This is the first time 
that clear linear behaviour, with the expected values of k, has been observed in hole-doped 
cuprates. 

The value of i7c2(0) ~ Hd 24 T obtained from these GL analyses may refer to a low 


field, unreconstructed Fermi surface. For fields greater than 24 T, we may be_o 
behaviour of weak superconductivity arising from the small electron pockets 


Dserving MF 


27| resulting 


from the appearance of CDW order. The GL type theory we applied assumes 5-wave su¬ 
perconductivity and we cannot rule out possible d-wave effects on the determination of iLc 2 - 
An obvious possibility is the Volovik effect whereby the Cooper pairs near the nodes on the 
Fermi surface are broken up, and consequently, A and k would increase. 

Alternatively, the diamagnetism that we observe above 24 T could be caused by super¬ 


conducting fluctuations. The Oil data in the insert to Fig. 3(a) show that it is ~ —100 A/m 
between 35 and 63 T. This is 5 times smaller and falls more quickly with field than predicted 
by theory 28(| for a 2D 5-wave superconductor at low temperatures and high fields. This is a 
robust statement because in the clean limit all parameters in the theoretical expression 281] 



















for Mc{B) above Hc 2 are known. Nernst data 29|] for OVIII crystal show satnration near 
30 T to the negative value expected for an electron pocket. This does not necessarily rule 
out bulk superconductivity above 30 T because in the presence for a CDW, the vortex core 
entropy - which dominates the Nernst effect - could be reduced. However at a qualitative 
level, the Nernst data between 24 and 30 T may be more consistent with superconducting 
fluctuations. Since torque magnetization is sensitive to superconducting fluctuations while 
thermal conductivity sees only the normal quasi-particles which are the only source of en¬ 
tropy, this may explain why we do not observe the sharp transition at Hk seen in Ref. 

0 . 

Finally, the diamagnetism of —90 A/m observed at 63 T might be caused by a reduction 
in ^in susceptibility associated with the stronger (3D) CDW order that sets in above 15 


ll| . The change required would be 1.36 A/m/T in Xd( 0). This is consistent with the 


signihcant decrease in diamagnetism between 40 and 50 K s 
the region where the 3D CDW seen at high helds goes away 


lown in the inset of Fig. 3(a) 

3 . 


The above discussion highlights the importance of competing CDW and superconduc¬ 
tivity instabilities 3, 30|. Little is known about the size of the CDW energy gap, or the 
MF behaviour expected for a d-wave superconductor just below Hc 2 as T —)■ 0 K. Therefore 
the linear H dependence of dM/dB we observe below might be a fundamental property 
of a d-wave superconductor. In other words, because of Volovik-type pair breaking effects, 
the MF transition at Hc 2 could have a discontinuity in d'^M/dB^, rather than in dM/dB, 
which is the standard MF result for the second order transition in a conventional s-wave 
superconductor. 

In summary, we observe diamagnetism in OVIII YBCO at helds up to 33 T and Oil 
YBCO at helds up to 65 T using torque magnetometry. The analysis uses a diherent model 
for the high temperature normal state susceptibility that gives a smaller correction at low 
temperature compared with earlier models. We also hnd that dM/dB departs from a linear 
lower held behaviour at helds Hd ~ idirr(O) ~ 24 T, and approaches a constant value at higher 
helds. We propose several candidate scenarios: a competing order scenario where a fully- 
hedged CDW at high held mostly suppresses the superconductivity so that the diamagnetism 
at high held could be attributed to bulk superconductivity; or a huctuation picture in which 
for H > Hd, the system crosses over to superconducting huctuation behaviour; and lastly, 
the apparent diamagnetism could rehect a reduction in spin susceptibility associated with 
the stronger 3D CDW order. It would be of interest to develop d-wave expressions for the 
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FIG. 5. (Color online) for both Oil and OVIII crystals show similar temperature dependences. 
Exponential fits to of 011(23.2 exp(—T/13.5))) and OVIII(23.7exp(—r/20.5)) give extrapo¬ 
lated values i7irr(0) = 23.2 and 23.7 T. These values are close to the low temperature for both 
crystals. Note that Hc 2 from GL fits (see main text) connects smoothly to 

MF magnetization and for the fluctuation contribution in the low temperature, high held 
regime, for comparison with our data. This could settle the debate over the existence of the 
high held vortex liquid region. 
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